Cosmic optical activity from an inhomogeneous Kalb-Ramond field 
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The effects of introducing a harmonic spatial inhomogeneity into the Kalb-Ramond field, inter- 
acting with the Maxwell field according to a 'string-inspired' proposal made in earlier work are 
investigated. We examine in particular the effects on the polarization of synchrotron radiation from 
cosmologically distant (i.e. of redshift greater than 2) galaxies, as well as the relation between the 
electric and magnetic components of the radiation field. The rotation of the polarization plane of 
linearly polarized radiation is seen to acquire an additional contribution proportional to the square 
of the frequency of the dual Kalb-Ramond axion wave, assuming that it is far smaller compared to 
f*) ■ the frequency of the radiation field. 
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p ■ I. INTRODUCTION 

The behaviour of electromagnetic waves in a curved background spacetime with torsion and its cosmological con- 
QO ■ sequences, has been an area of some interest in recent years. This is in view of its prospective implications for the 
low energy approximation to string theory. One way to investigate these implications is by identifying spacetime 
torsion yj, with the massless antisymmetric second rank tensor field, i.e., the Kalb-Ramond (KR) field, existing in 
most supergravity theories and as such in the massless sector of the most viable string theories . Certain physically 
observable phenomena result from the above analysis: a cosmic optical activity involving the rotation of the plane 
^y-j ■ of polaization of linearly polarized synchrotron rotation from high redshift galaxies investigated recently || being an 
example. The angle through which the polarization plane rotates is shown to be proportional, to leading order in 
, inverse conformal time (which decreases with redshift), to the rate of change of the KR dual axion field. It is also 
independent of the frequency of radiation. This latter aspect is in contrast to the well-known Faraday rotation, and 
, hence a new phenomenon. 

Another point to note is that the KR field, argued to be responsible for the effect, has been treated in as a 
perturbation on the Maxwell equations in a standard Friedmann-Robertson- Walker (FRW) cosmological background 
with both matter and radiation domination. Thus, it is assumed to have a negligible effect on shaping cosmological 
O i| background spacetime. One way of thinking about this is to imagine that the KR axion decouples from the radiation 
D . or matter (dust) fluid shaping cosmic geometry far prior to dust-photon decoupling, leaving behind a 'Cosmic KR 
1 Background' which affects incoming radiation from distant galaxies, albeit rather softly. As the universe expands 
further, this effect will no doubt gradually subside. The point made in || is that the effect may yet be observable in 
this epoch. This viewpoint has received support from [|| where the the effects of a time-dependent dilaton field are 
additionally incorporated, while demonstrating that the earlier findings have a degree of robustness. 

However, in these earlier assays, it is assumed that the axion and dilaton fields are spatially homogeneous, depending 
only on the conformal time coordinate rj. This is possibly quite justified considering the overall homogeneity of the 
universe over cosmological distance scales. 

In this paper, we generalize the scenario in Q by introducing a spatial inhomogeneity of a particular type into 
the KR axion field H: the massless Klein-Gordon equation obeyed by the axion quite naturally leads us to take the 
spatial dependence in the form of a plane wave propagating in space. However, taking into consideration the fact that 
the effects caused by the KR field should be confined to a very feeble disturbance on the overall homogeneity of space, 
we assume the the frequency of the axionic wave to be far smaller compared to that of the electromagnetic wave. 
This modification produces some interesting features, such as, it alters the mutual orthogonality of the electric and 
magnetic field vectors while inflicting a change on the Poynting equation as well. Furthermore, the inhomogeneity 
produces an additional rotation of the plane of polarization of radiation over that found in |3|]. This additional 
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contribution turns out to be porportional to the square of the frequency of the axion field while being independent of 
the wavelength of the radiation. 

The paper is organized as follows. For completeness and as a background, we briefly recapitulate in Section 2 the 
basic tenets of S, which lead to a modified set of Maxwell equations. The harmonic spatial dependence assumed for 
the axion field further modifies these equations as well as the one for the H field itself. Demanding that the waves 
retain their forms obtained in || in the limit H becomes space-independent, we obtain the equations representing 
circularly-polarized states. Carrying out a standard WKB type procedure, we solve these equations and calculate the 
rotation angle of the plane of polarization of the fields for a flat background spacetime in Section 3 and for a spatially 
flat spacetime in Section 4, considering separately the radiation and matter dominated cases therein. We make a few 
concluding remarks on our results in Section 5. 



II. GAUGE INVARIANT EINSTEIN-CARTAN-MAXWELL-KALB-RAMOND COUPLING 

A. Modified field equations 

The action for gauge-invariant Einstein-Cartan-Maxwell-Kalb-Ramond coupling is taken to be of the form 
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S = d x\/—g 



R(9,T) 1 1 ~ - x 

re 4 M -^W? 



(1) 



R{g, T) being the scalar curvature for the Einstein-Cartan spacetime where the connection contains the torsion tensor 
Tafiv (supposed to be antisymmetric in all its indices) in addition to the Christoffel term; re = 16ttG is the coupling 
constant; and H^x is the KR field strength three-tensor modified by U(l) Chern-Simons term arising form the 
quantum consistency of an underlying string theory: 



H^x = d^B^x] + -^VkA^F u x], 



(2) 



with, B v x being the antisymmetric KR potential which is considered to be the possible source of torsion. R is related 
to the scalar curvature R of purely Ricmannian (torsion-free) space-time by 



R{g, T) = R(g) + T^xT^ 



(3) 



The fact that the augmented KR field strength three tensor plays the role of spin angular momentum density (which 
is the source of torsion Q|) can be evidenced directly from Eq.(l) where the torsion tensor T^x, being an auxiliary 
field, obeys the constraint equation 



T, 



(4) 



Substituting Eq.(4) in the action (1) and varying the latter with respect to and respectively, two sets of field 
equations are obtained 



D^ uX = 



1 



-.d^V^gH^) = 



(5) 



and 



= ^H^ x Fx^. 



In addition, there is also the Maxwell - Bianchi identity 



Dfj, *F» V EE 
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(G) 



(7) 



Note here that all covariant derivatives are defined with the Christoffel connection and the Maxwell field strength is 
the standard 2-form F = dA. 

Now, expressing the KR field strength three tensor H^x = d\u.B v \] as the Hodge-dual to the derivative of the 
spinless pseudoscalar field H (the axion): 
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H^x = e% x D p H. (8) 

and substituting in Eqs.(6) and (7), the modified generally covariant Maxwell's equations are obtained in three- 
vectorial form 

D E = 2Vk Dff • B (9) 
£> E - D x B = [ D H B - DH x E ] + 0(k) (10) 

DB = (11) 
D a B - D x E = (12) 

where Z? M stands for the covariant derivative. On dropping all the higher order terms, as a first approximation, and 
retaining terms only of the order of ^/k, in a spatially flat isotropic FRW background with metric 

ds 2 =i? 2 (?7)(dr/ 2 -dx 2 ), (13) 

the set of equations (9)-(12) take the form 

V • E = 2VH • B (14) 

djt - V x B = -2[ d n HB -VffxE] (15) 

V-B = (16) 

d v B + V x E = (17) 

where r\ the conformal time coordinate, defined by dr\ = dt/R(t), R is the cosmological scale factor; and E = R 2 ~E 
and B = R 2 B. H is redefined by absorbing the y/R in it. 

It is easy to show from the very form of the KR field strength, viz., H^x — d^B^x], that it satisfies the Bianchi 
identity 

e^d p H^x = (18) 
which immediately implies that H satisfies the wave equation 

D p D p H = (19) 
In isotropic spatially flat universe this equation reduces to 

(9 2 - V 2 )H = 2^H (20) 

where H is taken to be a general function of both space and time coordinates, and the over-dot implies partial 
differentiation with respect to r/. Assuming a general wave solution for H, viz., 

H( V , X ) = H ( V ) cos p-x (21) 

Eq.(20) enables us to get an equation for H 

H + + p 2 H = (22) 

R 

We should point out here that as a consequence of our prior assumption that the overall homogeneity of the universe 
over long distance scales is not much disturbed by the inclusion of the spatial part in H, we are taking p to be much 
less compared to the wave number for the electromagnetic radiation. 

B. Modifications in electro-magnetic orthogonality and the Poynting equation 

From the field equations (14) - (17) we derive the wave equations for the electric and magnetic fields 

□ B = ( d 2 - V 2 ) B = 2 V x ( HB) - 2 V x ( VH x E ) (23) 

□ E = ( <9 2 - V 2 ) E = 2 V x ( HE) + 2Vi?xE - V(V-E) - 2 HB (24) 
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These equations indeed reduce to the pure Maxwell equations in the limit H — > 0, or a constant. Treating the 
axion H as a tiny purturbation over the Maxwell equations we argue that the solutions should have a form not much 
departing from the usual plane wave structure with the wave vector k perpendicular to both the electric and magnetic 
vectors. Moreover, the field equations (14) and (15) enable us to derive 

VH • B = (25) 

which, in view of the specific form of H, viz., Hq(ji) cos(p • x), implies 

p • B = (26) 

provided H ^ 0, which is the general case we are handling. 

This orthogonality of p and B [Eq.(26)] makes it easier to assume, for simplicity, that p can be taken to be 
orthogonal to E as well, i.e., p is either parallel or antiparallel to k. This is fairly justified, as it looks, from the 
similar wave nature of the electric and magnetic fields, at least in the limiting pure Maxwellian case as we are treating 
the modification caused by the KR field as a small purturbation over the Maxwellian behaviour. In fact, we may 
choose to observe the electromagnetic radiation which is travelling in the direction of propagation of the KR field. 

Considering the z-direction to be the propagation direction of the electromagnetic waves and as such for the axion 
(following the abovementioned assumption) we reduce the four-dimesional problem to a two-dimensional one with r\ 
and z being the only variables. The field equations now reduce to simpler forms 

□B = 2 HV x B + 2d v (H'e z x B) + 2 H"~E (27) 
□ E = 2 HV x E + 2d ri (H'e z x E) - 2 HB (28) 

where the over-dot and prime denote respectively the partial differentiations with respect to rj and z; and e z is the 
unit vector along the z-direction. Now, it is easy to show from the field equations (15) and (17) that 

d v (E • B) = B-VxB -E-VxE - 2 H B 2 + 2Vff-ExB (29) 

Considering the limiting plane wave behaviour of the solutions of the electromagnetic wave equations as H — ► , 
the magnetic and the electric fields can be expressed 

B(n,z) =B (r?, z) e- lkz (30) 
E(r,,z)=V (v,z)e- ikz (31) 

When the KR field actually vanishes, Eo, Bo = constant vectors x e zkv in the plane wave solutions of the pure 
Maxwell equations. Eq.(29) then gives E-VxE = B-VxB = 0, i.e., <9,,(E • B) = 0. Moreover, since B is 
in the direction of V x E, as is evident from Eq.(17), it follows that E • B = 0. When the KR field is present, but 
is only time-dependent (the case in ||), the vectors Eo and Bo in the solutions (30) and (31) again depend on time 
only and Eq.(29) reduces to 

d v (Eo-Bq) - - 2 H B 2 . (32) 

Clearly, E • B 7^ 0, which implies that the mutual orthogonality of E and B is lost. In the most general case 
where the KR field depends on both space and time coordinates, Eo and Bo are also spacetime-dependent and they 
satisfy the relation 

5„ (E • B ) = B • V x B - E • V x E - 2 H Bq + 2 VH ■ E x B (33) 

Here we find, one cannot ascertain conclusively that E ■ B is manifestly zero. It is, in fact, more justified to think 
that Eo ■ Bo is essentially non-zero in general, at least, by looking at the limiting behaviour when the KR field 
is stripped off the spatially depending part. The last term in Eq.(33), appearing due to the inclusion of the spatial 
dependence in H cannot, in general, compensate for the term which actually renders Eo • Bo non- vanishing in Eq.(32), 
as the temporal and spatial components of the KR field are completely separate entities. 

The Poynting equation in the present scenario can be obtained directly from the field equations (14) - (17). It is 
given by 

V • S + u em = - 2H E • B (34) 
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where S = (E x B) is the Poynting vector and ui em = |(E + B ) is the electromagnetic energy density. 
The distinction of this equation over the Poynting equation in the pure Maxwellian case is the presence of the term 
on the right hand side which is, in general, non-zero for reasons discussed above. In fact, this term vanishes in the 
limit H — > 0, i.e., when H becomes purely space-dependent. This is quite obvious since spatial inhomogeneity alone 
cannot bring in any change in the power conservation equation. 



C. Polarization States and Duality transformation 

Following H and rearranging terms of the components of the wave equations (27) and (28) (following the procedure 
in we obtain the following equations for the polarized states 

b± T 2iH'b± T 2iH'b± = b± ± 2iHb'± + 2H"e± (35) 

e± T 2iH'e± T 2iH'e± = e± ± 2iHe'± - 2Hb±, (36) 

where 

b±(v, z ) = B x (rj,z) ± iB y (r],z) 

e±(v, z ) = E x (i],z) ± iE y (r],z). (37) 

Note that the Eqs.(35) and (36) are converted into each other by the transformation e± — > b±, b± — * — e± provided 
the equation DH = H — H" = is obeyed. This is the usual electro-magnetic duality symmetry of the Maxwell 
equations. The Maxwell-KR system indeed possess this invariance in a flat spacetime background with cosmplogical 
scale factor R = 1, as is evident from Eq.(20). In a curved spacetime background, however, one does not have this 
invariance any more. 

Rewriting the Eqs.(35) and (36) as follows: 

b± T 2iH'b± T 2iH'b± = b'± ± 2iHV± + 2H"e ± = a±(<n,z) (38) 
e± T 2iH'e± T 2iH'e± = e± ± 2iHe' ± - 2Hb± = /3 ± (rj,z), (39) 

We seek the appropriate forms of the functions a(rj, z) and (3(rj, z) by examining the limiting forms of the above 
equations: 

I. Limit H' -> : 

In this limit, where H becomes purely a function of 77, the equations (37) and (38) reduce to 

b±= ft" ± 2iHb± = a±(H'->0) (40) 
e± = e± ± 2iHe'± - 2Hb± = 0±(H' 0) (41) 



Assuming solutions of the form 



we find 



b±(v,z)= b±( V ) e- lkz (42) 
e±(v,z)= 4M e- lkz (43) 



b± = - (k 2 T 2kH)b± = a±{H' -> 0) (44) 
e'±= — {k 2 =F 2kH)e± - 2Hb± = P±{H' -> 0). (45) 

II. Limit H -> : 

This limiting case implies H to be a function of z only but as is evident from Eq.(22) H' is merely a constant. 
Therefore Eqs.(37) and (38) reduce to 

b± T 2iH'b±= b± = a±(H^0) (46) 
e± T 2iH'e±= e± = P±{H -»• 0) (47) 

Again assuming solutions of the form 
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we get 



M»7>*) = bt(z)e ik *> 
e ± ( V ,z)= ef(z) e lfe " 

b'l= - (fc 2 T2kH')b ± = a±(H^Q) 
4 = - (fc 2 T2kH')e± = f3 ± (H -> 0). 



(48) 
(49) 

(50) 
(51) 



By looking at these limiting forms of a± and /3± given in Eqs.(44),(45) and in Eqs.(50),(51) it seems reasonable to 
suggest the following simplest possible structures of these functions: 



a±(ri,z) 



fc 2 =p 2k(H + H') b±(ri,z) 

k 2 =p 2k(H + H')] e±(ri,z) - 2Hb±(r),z). 



Moreover, setting 

e±(r),z) = a±(rj,z) b±(rj,z) 
we write the equations (37) and (38) in a more elegant form 

b± T 2zH'b± + \k 2 =f 2k(H + H') t 2iH']b±(r),z 



and 



bl ± 2iHb', 



e± T 2iH e± 



k 2 T 2k(H + H') + 2H"a± b±(r],z 



k 2 =f 2k{H + H') t 2iH' + 2 



H 

a± 



e±{n,z) = 



4 ± 2iije' ± 



fc 2 =f 2k(H + H') e±(r),z) = 



(52) 

(53) 
(54) 

(55) 

(56) 
(57) 

(58) 
(59) 



III. FLAT SPACETIME BACKGROUND 

In order to get a preliminary idea as to how the coupling of a spacetime dependent KR field to Einstein-Maxwell 
theory affects the electromagnetic waves, and thereby effects in an optical activity in the radiation coming from 
distant galactic sources, we consider the simplest situation — that is, of a flat universe with cosmological scale factor 
R(rj) = 1. Admittedly, quantitative details of results of this section are cosmologically untenable for obvious reasons. 

The equation of motion (22) for H can be solved readily to obtain 



H(rj,z) = (cisiripr] + C2COsprj) cospz 



(60) 



where ci and C2 are arbitrary integration constants. Demanding that the above solution must reduce to the form 
(hrj + ho) which is the solution of Eq.(22) in the limit p — > we infer c\ — h/p and C2 = ho. Here h and ho are 
the same arbitrary constants denoted in B. Making a Taylor series expansion of the various functions appearing in 
Eq.(60) around p = we write 



ff(7?,z) = {hrj + 



oy 2 V 3 



+ M + hr)Z z + h a z z + 0(p 4 ) 



Substituting this H in Eqs.(55) and (57) and assuming solution of the standard WKB type given by 

b±(v,z) = b e lk s ±^ z > 



(61) 



(62) 



with 
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S. 



2 

k 2 



we get after partial integrations of Eqs.(56) and (57) with respect to r/ and z respectively 



b±(i]i z ) = b explikr/ =p i 



p 2 / /l7y 

^ - y (— + M 



hz r\ 



and 



b±{rj,z) = b exp < ikg±(r]) — ikz =F i 



P -{h V + h )z 2 ) + 0(p 4 ) 



ikf±(z) 



(63) 



(64) 



(65) 



Comparing these two expressions we assert the forms of the arbitrary functions f±{z) and g±{ii) and write b±(r],z) 
as follows: 



b±(i]-, z ) ~ b exp i.ik(rj — z) =p i 
A similar approach for e±(j], z) yields 

e ±( f l; z ) — 6 exp < i/c(?7 — z) =p i 



/;r/ 



/;?/ 



p 2 / hrj 



+ h Q rf - h z z + 0(p 4 ) 



2 



/i?7 3 
~3~ 



h T] 2 



0{ P A 



o 



o 



(66) 



(67) 



which differs from Eq.(66) only in the constant coefficient e and in the higher order O (■£)• However, it should be 
mentioned here that while using the WKB technique we are assuming that the function a± is not increasing rapidly 
as k increases. In fact, in absence of the KR field, when we have the plane wave solutions of Maxwell's equations, 
a + = a_ = e/b = constant. The spacetime dependence of a± comes only in presence of a spacetime-dependent 
H. Since H is being treated as a small purturbation over the Maxwell field, it is rather plausible to think a± to be 
not much different from the constant e/b and a very slowly-varying function of spacetime. But the constant e/b is 
arbitrary and cannot generically be argued as increasing with k. Therefore, the assumption that a± remains quite 
invariant as k increases is fairly justified. Using WKB method, the solutions obtained above involves a± only in the 
higher order O (£). 

Now, the circular polarization states are defined by b± and e± and the extent of the optical birefringence due the 
presence of the KR field can be estimated directly by calculating the rotation angle of the plane of polarization of the 
electromagnetic wave, which is given by the phase difference (f> m ag = \[ ar 9 b+ — arg 6_] for the magnetic field and 
(j> e iec = \[ ar 9 e + ~ ar 9 e_] for the electric field. The phase shift is given by 



<f>mag(W, %) ~ 4>elec(r), z) 



p 2 ( hvf^ 
hv + y ( -g 1- h v 2 



h z' 



(68) 



for h,p << k. 

It is interesting to see that the change in the rotation angle, calculated here, over that found in || for flat universe, 
is primarily given by the p 2 -dependent part. But as p is considered to be very small we infer that this change is rather 
insignificant. 

IV. SPATIALLY FLAT FRW SPACETIME BACKGROUND 

We now turn to less trivial background spacetimes. We consider a spatially flat expanding universe dominated by 
radiation and matter, in turn. 



A. Radiation dominated Universe 



The scale factor for this model, in real time, is given by 

[R(v)} RD 



v_ 

rj r 



(69) 
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^ In 

where i] r = (87rGeo/3) , eo being the primordial radiant energy density. 
Substituting this in the equation of motion (22) for H we obtain 



jfH + 2r]H + p 2 rj 2 H 







which has the form of a transformed Bessel equation with solution 

H(r],z) = 77 * [ci J 1/2 (pv) + c 2 Y 1/2 (pr])] . 
Simplifying the Bessel functions of first and second kinds, viz., J and Y, the above solution can be written as 



(70) 
(71) 



H(rj, z) = — (ci sin pi] + c 2 cos prj) cos pz. 
Imposing again the boundary condition that this must reduce to the limiting form 



ofEq.(22) — as p — > , we set c\ = y- and c 2 = — h r , with h r = h r/ 2 , whence 

H(n,z) = — (— sin pri — h r cos pn ) cos pz. 
V V P J 

Plugging in the Taylor expanded form of this H in Eqs.(56) - (59) and using the same WKB technique as for the 
flat universe, we obtain 



(72) 

the solution 
(73) 



b±(rj, z) — b exp 



ik(r] — z) ± i 



h r 

n 



P 2 ( hoV 2 



h r r] - h z 2 + 0(p A ) 



and similar expression for e±(i], z). 

The phase shift in this case is given by 

(T),Z) W (f>elec(V,z) 

for h,p « k and h r = h r] 2 . 



K P f, h rf 2 



(74) 



(75) 



B. Matter dominated Universe 

In this case, where the scale factor is governed by 

[R(v)] MD ^ 



with rj m = (8nGp /3) (p - the initial matter density), Eq.(22) reduces again to a transformed Bessel equation 



(76) 



rj*H + 4r)H + p 2 r) 2 H = 



having simplified solution 



1 



H(r], z) = — [(ci - c 2 prj) sin pr\ — (c 2 + c\ prf) cos prj] cos pz. 



n 



(77) 



(78) 



Determining the constants c\ and c 2 using, as before, the boundary condition on H in the limit p — > , we write 

K 

where h m — hrf. With this H, we obtain using the WKB method 



H(r], z) = — (cos prj + pr\ sin prj) cos pz. 



(79) 



b±(i], z) = b exp 



ik(r] — z) ± i 



and similar expression for e±(r), z). The phase shift can be calculated 



3t7 c> \r/ 6 



(80) 



(81) 



for << k. 
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V. CONCLUSIONS 



One rather surprising aspect of our finding is the loss of orthogonality of the electric and magnetic vectors in the 
radiation field which exists even in the case of a spatially homogeneous axion field, so long as the axion is time- 
dependent. This will indeed affect the measurement of the precise rotation of the plane of polarization, although for 
large redshift sources to which we confine, this effect may be ignored for all practical purposes. 

In general the electric and magnetic vectors seem to have solutions generically represented as 

(Electric/ Magnetic field combinations) = (constant) e l(kz ~ uv) e lH{v < z) (82) 

The additional contribution to the angle of rotation arising out of the spatial inhomogeneity introduced in this paper 
is actually quite small under our assumption that the wavelength of the axion wave is far larger than the radiation. 
Thus from an observational standpoint this new effect is not too significant, although it is still quite distinct in its 
behaviour from Faraday rotation. 

We have used standard WKB type methods to arrive at some solution to the complicated set of equations which 
arise even at the lowest order. Among other things we have explored the properties of the electric and magnetic fields 
of the Maxwell-KR system. We also checked the validity of the standard electric-magnetic duality and pointed out 
the existence of a very simple solution to the complicated equations. 
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